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An inertia-gravity wave (IGW) propagating in a vertically sheared, rotating stratified 
fluid interacts with the pair of inertial levels that surround the critical level. An exact 
expression for the form of the IGW is derived here in the case of a linear shear and 
used to examine this interaction in detail. This expression recovers the classical values 
of the transmission and reflection coefficients |r| = e~'^^ and \R\ = 0, where = 
J(1 -I- — 1/4, J is the Richardson number and v the ratio between the horizontal 

transverse and along-shear wavenumbers. 

For large J, a WKB analysis provides an interpretation of this result in term of tun¬ 
nelling: an IGW incident to the lower inertial level becomes evanescent between the 
inertial levels, returning to an oscillatory behaviour above the upper inertial level. The 
amplitude of the transmitted wave is directly related to the decay of the evanescent 
solution between the inertial levels. In the immediate vicinity of the critical level, the 
evanescent IGW is well represented by the quasi-geostrophic approximation, so that the 
process can be interpreted as resulting from the coupling between balanced and unbal¬ 
anced motion. 

The exact and WKB solutions describe the so-called valve effect, a dependence of the 
behaviour in the region between the inertial levels on the direction of wave propagation. 
For J < I this is shown to lead to an amplification of the wave between the inertial levels. 
Since the flow is inertially unstable for J < 1, this establishes a correspondence between 
the inertial-level interaction and the condition for inertial instability. 


1. Introduction 

Inertia-gravity waves (IGWs), that is, internal gravity waves with frequencies close 
enough to the Coriolis frequency to be affected by rotation, are ubiquitous in both the 
ocean and the atmosphere (for recent observations see Marshall et al. 2009 and Hertzog 
et al. 2008, respectively). In the ocean, they are efficiently excited by surface winds es¬ 
pecially in the region of the atmospheric storm tracks (Alford 2003). In the atmosphere, 
external sources such as mountains and convection produce IGWs (Scavuzzo et al. 1998, 
Lott 2003); internal sources, like those appearing during the life cycle of baroclinic in¬ 
stabilities are also efficient (Plougonven and Snyder 2007, Sato and Yoshiki 2008). The 
impact of atmospheric IGWs on the global climate is now well established (Andrews et 
al. 1987), and a current challenge is the quantification of the non-convective sources that 
are parameterized in middle atmosphere models (Zuelicke and Peters 2008, Richter et 
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al. 2010). In the ocean, IGWs are important for several processes including small-scale 
mixing and energy dissipation (Ferrari and Wunsch 2009). 

IGWs are generated in or propagate through regions where intense jets exist (Whitt 
and Thomas 2013), and they interact with the jets. These interactions involve critical 
levels, that is, regions where the Doppler-shifted frequency of the wave is zero (Booker 
and Bretherton 1967) or, when rotation is taken into account, inertial levels where the 
Doppler-shifted frequency equals the Coriolis frequency (Jones 1967). Such interactions 
arise, for instance, when a low-level front passes over a mountain ridge, yielding direc¬ 
tional critical levels almost everywhere at low altitudes (Shutts 2003, Shen and Lin 1999). 

A central result, valid both with and without rotation, is that an IGW originating from 
z ^ —oo, propagating though critical and inertial levels and radiating away as 2 —>■ cxd, 
is not reflected and is absorbed by a factor 

|r| = exp(—TT^), where ^ = \/j(l -I- v'^) — 1/4, (1.1) 

J > 1/4 is the Richardson number, and v = l/k is the ratio of the horizontal components 
of the wavevector respectively transverse to and aligned with the shear. This result is 
obtained for an inviscid, hydrostatic fluid under the assumption of uniform vertical shear 
and Brunt-Vaisala frequency. It gives the ratio of the amplitude of the transmitted wave 
to the incident wave and, in the rotating case, conceals a complex behaviour between 
the inertial levels. This is exemplified by the fact that the change of the wave amplitude 
across the lower and upper inertial level are given by factors 

exp (ttz/) and exp (— ttz/) , (1-2) 

respectively (see Grimshaw 1975), and by the fact that the wave amplitude does not 
change at the critical level. While they cancel overall, the factors in & reveal a strong 
dependence of the wave amplitude between the two inertial levels on the direction v = l/k 
of the wave vector. This phenomenon, referred to as the “valve effect”, has been con¬ 
sidered in earlier papers (Grimshaw 1975, Yamanaka and Tanaka 1984), but none are 
entirely satisfactory since they do not provide exact solutions or approximate solutions 
over the entire vertical axis. In this respect it is worth recalling that the overall transmis¬ 
sion in cu is the same with or without rotation (Jones 1967) because the solutions as 
z —>■ ±c» can be connected by integrating the relevant differential equation along a path 
with \z\ » 1 on which the solution is asymptotically unaffected by rotation (see Fig.[T]). 
This is made possible by the structure of the Stokes lines which, unlike in the familiar 
Airy equation for instance, do not go to c», but simply join the two inertial levels along 
the real z-axis. One of our aims is to give an exact solution for z = 0(1) that displays 
a dependence on rotation and zz; another is to provide a description of the solution that 
reconcile the absorptive properties in CH) and m- 

Another motivation for the paper is the relation between critical and inertial levels, 
and instability. This relation is well known: classical shear instabilities, for instance, can 
produce internal gravity waves (Lott et al. 1992, Lott 1997), in which case a critical 
level is necessarily present (Miles 1961); conversely, the interaction at a distance between 
gravity waves helps explain stratified shear flow instability (Rabinovitch et al. 2011). 
Similarly, recent works on the coupling between balanced waves and unbalanced waves 
show that critical and inertial levels are associated with unbalanced instabilities (Sakai 
1989, Molemaker et al. 2005, Plougonven et al. 2005, Vanneste and Yavneh 2007, Sutyrin 
2008, Gula and Zeitlin 2010). In the ocean, these unbalanced instabilities may have an 
important role in providing dynamical routes to dissipation at small scales (McWilliams 
2003). 

In this context of flow stability, the fact that (HH) remains valid in the rotating case 
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is surprising. Indeed, in the non-rotating case the conditions for gravity-wave absorption 
and flow stability are the same, J > 1/4 (Miles 1961, Howard 1961), leading to an in¬ 
terpretation of instability in terms of gravity-wave amplification (Lindzen 1988). In the 
presence of rotation, however, symmetric instability (with k = 0) occurs for J < 1 (Stone 
1966, Bennetts and Hoskins 1979), but the critical value J = 1 does not directly relate to 
the wave-absorption properties in (HU or (HU. For J ^ 1 and with horizontal bound¬ 
aries, non-symmetric disturbances can compete with the symmetric ones (Mamatsashvili 
et al. 2010), but again the relation with the absorptive properties of the critical layer is 
not clear. One hint of a relation for symmetric and near-symmetric disturbances is the 
behaviour of the amplification factor in (11.21) ds v = l/k ^ oo, but this does not involve 
J. 

The purpose of this paper is to reconcile the absorptive properties in HU and (HU by 
giving an exact solution of the Taylor-Goldstein equation, with constant shear, stratifica¬ 
tion and without boundaries, that is valid over the entire range of altitude in the rotating 
case. The interpretation of this solution is facilitated by an asymptotic WKB treatment 
valid when J ^ 1. Our method is related to that of Lott et al (2010, 2012 hereinafter 
LPV12) who showed the importance of inertial levels for the spontaneous emission of 
gravity waves by potential-vorticity (PV) anomalies. Specihcally, they showed that a PV 
anomaly in a shear induces a balanced response that decays exponentially in the vertical 
up to the inertial levels, then takes the form of a propagating IGW. 

In the present paper, we consider an IGW that is incident from below to a pair of 
inertial levels. It changes nature across the lower inertial level, taking an exponentially 
decreasing, balanced form, and is restored to an IGW of much weaker amplitude across 
the upper inertial level. Qualitatively, the amplitude of the transmission coefficient in 
HU is set by the value of the balanced solution immediately below the upper inertial 
level. Our exact and approximate results provide the properties of the solution for a 
broad range of values of J and give a complete description of the solution between the 
inertial levels, where the valve effect is manifest. For J < 1, in particular, we identify a 
form of disturbance amplihcation, whereby the Eliassen-Palm (EP) flux of the wave is 
amplified as the wave crosses the lower inertial level. This establishes a clear correspon¬ 
dence between the absorptive properties of the shear layer and the criterion J < 1 for 
flow instability. 

The plan of the paper is as follows. The exact solution for the vertical structure of 
the wave is derived in section 2. As the mathematics are quite involved but present only 
few conceptual difficulties much of the derivation is relegated to an Appendix. Section 
3 develops an interpretation of the exact result in terms of tunnelling using a WKB 
analysis for J ^ 1. The inertially unstable case J < 1 is analysed in section 4. The paper 
concludes with a brief discussion in section 5. 


2. Exact solution 

We start from the linearized hydrostatic-Boussinesq equations for the propagation of 
a three-dimensional disturbance in a uniformly stratified sheared flow Uq = (Az,0,0), 
where the shear A > 0 and the Brunt-Vaisala frequency N are constant. When the Cori¬ 
olis parameter / is also constant, a stationary monochromatic disturbance with vertical 
velocity of the form 



( 2 . 1 ) 
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with Lo the absolute frequency, satisfies the vertical structure equation 



(1 -f iy^)J 2iv 

U' e) ^ 

[ e ^ CM 


W = 0, 


where 


N'^ I kAz Lo 

e = —-y. 


( 2 . 2 ) 


(2.3) 


This is Eq. (6) in Yamaha and Tanaka (1984), but see also Inverarity and Shutts (2000). 
Note that in the absence of horizontal boundaries, as assumed here, and for k ^{), there 
are no unstable modes so we can take Imw = 0; for fc = 0, inertial instability occurs for 
J < 1, leading to growing modes with Imw ^ 0 which we discuss further in section 4. 

Equation (12.2|) has two singularities at the inertial levels ^ = ±1 that require a special 
treatment; the critical level ^ = 0 is an apparent singularity. Following LPV12, we write 
the exact solutions of (j2.2l) as 

lEi™)(e) forC>l, (2.4a) 


AWr\0 + for - 1 < ^ < 1, 

CWf\0 + DWi^\0 for^<-l. 


(2.4&) 
(2.4c) 

In (|2.4al) - (l2.4eD A, B, C and D are constants, and the W functions can be expressed 
in terms of hypergeometric functions (see Appendix, Yamanaka and Tanaka 1984 or 
Shutts 2001). In (|2.4al) we have retained the solution with the asymptotic form 

as ^ oo, (2-5) 

corresponding to an upward-propagating gravity wave (Booker and Bretherton, 1967). 
In (|2.4cl) the two functions are such that the asymptotic form is given by 

VE(i) asC^-oo. (2.6) 

To evaluate A, B, C, and D, we first use asymptotic expansions of (I2.4al) - (|2.46l) and 
(|2.46l) - (l2.4el) in the vicinity of ^ = 1 and ^ = — 1, respectively (see Appendix). Next we 
follow Booker and Bretherton (1967) and introduce an infinitesimaly small linear damping 
to determine the physically relevant branches of multivalued functions (see Fig. [1]). This 
is equivalent to shifting the path of integration of (j2.2ll into the lower half of the complex 
plane so that 

^-l = (l-^)e-“ for ^<1, and , . 

^ + 1 = + for ^<-1. ^ ’ 

Using this, the asymptotic expansions of (|2.4a|) and (I2.46|) can be continued for ^ 1“ 

and ^ — 1“ and compared with the corresponding expansions of (j2.46l) and (I2.4cp (see 
the Appendix for an illustration of the procedure in the simpler context of the WKB 
solution). Matching expansions near ^ = 1 gives 


a' = e-‘'^ (aA + a"B ), 

/?' = PA + P"B, 


( 2 . 8 ) 


where the coefficients a, /3, a',etc are known functions of /r and v given in the Appendix. 
Similarly, matching near — 1 gives 


I3A-P''B = {P"'C + I3'D), 


(2.9) 


a'"C + a'D = aA- a''B. 

The constants A, B, C and D are then deduced by solving the linear system dini-diii). 
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infinitesimal paths WKB paths ^ f 




-oo 


Figure 1. Paths used for the integration of the vertical structure equation (1^ : see also 
Figure 1 in Jones (1967) and note that our real ^-axis is vertical while Jones’s Is horizontal. 


Note that the valve effect, involving an exponential dependence on v, can be traced 
to the factors in (j2.8l) - (j2.9ll . with the amplification across the lower inertial level 
associated with the factor in ()2.9I1 and the balancing attenuation across the upper critical 
level associated with that in (12.81) . 

After involved manipulations it can be shown that the reflection and transmission 
coefficients are given by 

i?=^=0 and T=^=ie-f^^, (2.10) 

where = J(1 + z/^) — 1/4 (see Appendix A.l). This is the result of Jones (1967), which 
can be obtained much more directly by integration along a semi-circle with |^| ;g> 1, as 
recalled in the introduction (see Fig. [1}. Our solution provides the details of the solution 
between the inertial levels that this result ignores. In particular, it makes it possible to 
compute the EP flux, proportional to 

= Re j > (2-11) 


and constant apart from jumps at the inertial levels (this conservation law is derived 
directly from (12.21) by multiplying by W* and integrating by parts). We can obtain the 
three distinct values of the EP flux on either side of and between the inertial levels using 
the first terms in the expansion of the hypergeometric functions for ^ > ±oo and ^ 0. 

Computations detailed in the Appendix show that the ratios of the EP flux between and 
above the inertial levels to the EP flux below are 


^z(II) g27r!v _|_ ^ 

|^Z(I)| ~ Q^TTfJ. _ 


^z(III) 


and 


-2TTfl 


( 2 . 12 ) 
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EP flux ratio: Log,o(F"'"VlF'''^l) 



Richardson Number (j) 


Figure 2. Ratio of the EP flux between the inertial levels to the incident EP flux as a 
function of the Richardson number J and wave aspect ratio u — l/k. 


where the absolute value in the denominators is required because the incoming EP flux 
< 0 whereas the other two fluxes are positive (see Appendix). The first of these 
ratios describes the crossing of the lower inertial level by the incident waves and displays 
the strong dependence on the direction u of the wavevector that characterises the valve 
effect. It is shown as a function of u and J in Fig. [51 The figure indicates that, for 
J < 1, the ratio can exceed 1 and can take very large values for u > 0. We discuss this 
phenomenon further in section 4. 


3. WKB approximation and tunnelling 


A WKB solution of (I2.2|l . strictly valid for /i ^ 1 but qualitatively correct for ^ ^ 1, 
sheds light on the behaviour of the solution and on the valve effect in particular. To 
derive this solution, it is convenient to follow Miyahara (1981) and recast (12.211 in the 
canonical form 


where 


4/55 + 


+ 1/4 
^ 2-1 


(3 + u2)^2_2- 
^2(^2 -1)2 _ 


4 - = 0 , 


W(0 = ?(? - 1)-1/2+W2 ^)-l/2-W2 ^)^)^ 


(3.1a) 


(3.1&) 


Assuming /r ^ 1, we then introduce the expansion 


and find at orders /i and 1 that 


(j) = —r^=, ^-0 = (^^ - 1)^^^ and 
VC -1 


.1 +(1/2/2-7/8)52 

ae - 1)'/" 


(3.2) 


(3.3) 


where e = ±1 selects the two possible solutions. The expansion (13.21) breaks down in 
the regions ||5| ~ 1| = surrounding the inertial levels where an approximation 

in terms of Hankel functions can be constructed (see LPV12 for the rescaling leading to 
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this conclusion). Here we avoid this complication by integrating (I3.1al) along a path that 
remains at a distance larger than from the inertial levels and avoids the branch cut 
(see Fig.[T]). This path does not cross the only Stokes line (where Ivaf = 0, e.g., Ablowitz 
& Fokas 1997), namely ^ G [—1,1]. Therefore, a single WKB solution, determined by the 
radiation condition to correspond to e = 1, can be continued from >ooto^—oo. 

To make the absorptive properties more transparent, we express the approximations at 
order /r for a unit-amplitude incident wave (rather than for a unit-amplitude transmitted 
wave in the exact case, see Appendix A.3 for details). This gives 

^ZMln(«+v^) (3.4a) 

g-/x(7r/2+sin-iO for-l<5<l, (3.46) 

fore<-l. (3.4c) 

Here, the function g{^) = -|- groups the factors unaffected by 

the absorptive properties of the shear layer; because it is expressed in terms of absolute 
values of 1, the complex powers —l/A±ii'/2 apply to purely positive real arguments, 
and the effect of the continuations appears explicitly through exponential factors in 

(I3.4al) -- (I3.46I) . Note that for |^| 1, k. and g{Ifj ~ sign(^)|^|^/^, 

so recovering the incident wave in (12.611 (up to an irrelevant phase 

factor). 

The WKB solution (13.41) is instructive in several respects. First, as ^ ±cxd, it repro¬ 
duces exactly the ratio between the amplitudes of the upward-propagating IGWs (12. 5|) 
and (12.611 found from the exact solution. In particular, the reflection and transmission co¬ 
efficients are the exact ones in (j2.10ll . The WKB solution also shows that the amplitude 
of the solution is dominated by the exponential term e“^®™ ^ in (13.461) that rapidly 

decays monotonically between the inertial levels. This decay is strongly reminiscent of 
the exponential decay that characterizes the neutral solution in the quasi-geostrophic 
approximation. This can be made more transparent by noting that near the critical level 
^ = 0, the 0(/i“^) term (retaining 4/]^ in (13.21) 1 is proportional to the quasi-geostrophic 
solution (1 -b (see Eq. (2.16) in LPV12). Away from the immediate vicinity of 

^ = 0, the decay is in fact faster than predicted by the quasi-geostrophic approximation 
since | sin“^ The valve effect is also evident in the WKB solution, through the 

factor in (13.461) . The dependence of the amplitude on n associated to this factor 

is however weak compared to the dependence on J and |i/| stemming from the factor 
g-M(7r/2+sin-io _ g-«) since J » 1. This factor represents a 

very rapid amplitude decay of the solution above ^ = — 1, and corresponds to an absorp¬ 
tion that increases rapidly with J and |jz|: the asymmetry that favours perturbations with 
z/ > 0 over those with z/ < 0 between the inertial levels is of little signihcance compared 
to this absorption. 

The WKB solution represents the exact solution accurately even for moderately large 
J. For J = 5, for instance, the two solutions almost coincide everywhere except in regions 
close to the inertial levels (see Fig. |3|). The WKB approximation deteriorates around the 
lower inertial level for |^ + 1| ^ 0.05 (see Fig. |31 d), corresponding to an 0(/i“^) distance 
to the inertial level, as expected. An analogous breakdown of the WKB approximation 
occurs for ~ 1, but in this case the solutions are very small. 

Figure [3] also illustrates how the solution behaves as an upward-propagating wave - 
with real and imaginary parts in quadrature - outside the inertial levels (see Fig. [3^), 


<KB = e- 

^ 5(0 

<1<B = 9(0 
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a)J=5,v=0.1, global view b)J=5, V=0.1, lowest inertial level 




Figure 3. Vertical velocity W for the exact solution (|2.4al) - (l2.4c|l (ReVF, black solid; ImW, 
black dashed) and the WKB approximation (I3.4al) - (l3.4el) (ReWwKs grey solid and ImVFwKB 
grey dashed) for J — 5, v — 0.1: (a) global view; (b) zoom on the lowest inertial level. 


and as an evanescent perturbation - with real and imaginary parts in phase - between 
the inertial levels (see Fig.jSJo). This behaviour can be interpreted as a form of tunnelling, 
analogous to, but different from the classical tunnelling between turning points as de¬ 
scribed, for instance, by Bender and Orszag (1978). The continuity of the arguments of 
the last exponential factors in (13.41) shows that the amplitude decrease of the transmit¬ 
ted signal by the factor e~^'^ can be immediately related to the decay of the solution 
W^wkb between the inertial level as ^ since this decreases by the same factor 

In the tunnelling interpretation, the valve effect seems to play no role, because the 
amplitude jump through the lower inertial level is exactly balanced by the inverse 

jump at the upper inertial level. 

An obvious difference with classical tunnelling is the absence of a fully conserved 
flux. In classical tunnelling, wave-action flux conservation leads to the relation \T\'^ + 
|i?^| = I that constrains the transmission and refection coefficients. Here, the analogous 
conservation is that of the EP flux, but because this can jump at the inertial level 
(Eliassen and Palm 1961, Grimshaw 1975), there is no associated constraint on T and 
R, making the absence of reflection i? = 0 compatible with partial (indeed very weak) 
transmission |r| < 1. Regarding the EP flux, we note that the WKB solution between the 
inertial levels (|3.4&P suggests a zero flux, since it is evanescent (see for instance Gill 1982). 
The flux is non zero, however, and can be captured by extending the WKB analysis to 
take into account the Stokes phenomenon that arises at ^ = 1. The relevant computation, 
carried out in LPV12, reveals the presence of an exponentially small solution increasing 
exponentially with ^ to be added to (I3.4&I) . Taking this into account gives a non-zero 
flux, deduced from Eq. (2.31) in LPV12 to satisfy 


^ 2 ( 11 ) 

consistent with (12.121) for ^ ^ 1. 




(3.5) 


4. Valve effect and inertial instability 

To describe further the valve effect, we return to the exact solutions and first present 
vertical profiles of W when J = 5 and for i/ = ±1 in Fig. 0]. The global views in Figs.|4^ 
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I ' I ' I '' III ' I I II ' I '' I \ li ' I 

-1 - 0,5 0 0,5 1 -1 - 0,5 0 0,5 1 

WandF' WandF* 


Figure 4. Vertical velocity W for the exact solution (I2.4al) - (l2.4cl) {\W\ thick grey solid, KeW 
grey solid, and ImW grey dashed) and EP flux from (12.1111 (black solid) for J = 5 and 
V = ±1. The EP flux is normalized by its incident value: (a) and (c) global views for v — ±1; 
(b) and (d) zooms on the lowest inertial level. 

and[4l: show that the solutions are almost identical. Even though a closer examination of 
the solution at the lowest inertial levels (Figs. |4 Jd and |4jl) shows that the solution with 
u > 0 is amplified whereas that with u < 0 is attenuated, this difference matters little 
because of the very rapid decay with ^ of both solutions between the inertial levels. Also 
in Fig.m we see that the EP flux is very small between the inertial levels consistent with 
expression (12.121) for the ratio of the EP flux between the inertial level to 

the incident EP flux. For J ;g> 1, this ratio is approximated by (ixa which shows, as 
discussed, that the valve effect plays only a minor role in modulating it. 

For J < 1, the situation is different, and the valve effect can lead to an intensification 
of the perturbation between the inertial levels. This is illustrated in Figure [5] which shows 
profiles of W and F^ for J = 0.5 and v = ±2. When u = 2 in Figs.[S^-b, the incident wave 
is clearly amplified as it passes through the lower inertial level (Fig.[5]3), and the solution 
between the inertial levels does not display the almost exponential decay with altitude 
that characterizes the solution when J ^ 1. The behaviour for v = —2 illustrated in 
Figs.Et-d is then completely different since the incident wave is almost entirely absorbed 
at the lowest inertial level. This pronounced difference between positive and negative v 
translates in the EP fluxes ratios in ()2.12D as Fig. [5] confirms. More importantly, (12.121) 
tells that the disturbance with v > \/{J — 1/4)/(I — J) are amplified between the inertial 
levels when J < 1 (in this case the ratio is larger than 1). More than this, for |u| ^ 1, 

pz(II) 

_ _ ^ p,2Tr{i^-Vj\v\) 

|F^(i)| 

so the amplification is arbitrarily large for u —>■ +oo when J < 1. 


(4.1) 
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Figure 5. Same as Fig. |4]but for J = 0.5, v = ±2. In a) and b), the EP flux is normalized by 
its amplitude between the inertial layers. Note the change of scale in panel (d). 

It is interesting to relate this result to inertial instability, which also occurs when 
J < 1 (Stone 1966), a condition that implies that the background flow potential vorticity 
is negative (Bennetts and Hoskins 1979). Inertial instability - or symmetric instability 
- arises for perturbations with fc = 0, that is, v = oo, a case that we have not treated 
explicitly. It is examined in Stone (1966) and can be recovered easily by replacing the 
non-dimensional variable ^ defined in (j2.3l) and clearly inappropriate for fc = 0 by the 
original z. This reduces the equation for W to 

{f - uj^) 114^ -h 2fA/1H4 - JIS^l^W = 0 (4.2) 

which admits the simple solution W oc exp(imz) for 

(jj2 2d -h 1) , (4.3) 

where 5 = (A/)/(/m). A purely imaginary frequency w, corresponding to symmetric 
instability, occurs for J < 1, with maximum growth rate Imw = f ^Jl/J — 1 for 5 = 
-IIJ. 

Although our analysis is not adapted to treat this case, because when fc —0 the 
inertial levels go to infinity and it becomes impossible to impose radiation conditions 
beyond them, the condition J < 1 for IGW amplification between the inertial levels 
establishes a clear connection between the valve effect and the inertial instability criteria. 
This connection with inertial instability is not the only one. For instance, the fact that 
waves with large v are amplified between the inertial levels suggests that the shear layer 
response to a random incident IGW field can become very large between the inertial 
levels and preferentially symmetric as are inertial instabilities. Also, the response in this 
transverse plane will be tilted in the direction of the isentropes (see Figs. 2d and 2h in 
LPV12) as is also the case for inertial instability since d < 0 for them (see (14.311 b Note also 
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that a stabilization of the flow by the disturbance is expected to occur in the nonlinear 
regime. Indeed the large increase of the EP flux at the lower inertial level accelerates the 
mean flow there, whereas the decrease at the upper inertial level decelerates it there. As 
a result, the positive shear tends to be reduced: as is often the case, the development of 
the instability weakens its cause. 


5. Conclusion 

The problem of the absorptive properties of gravity waves by critical levels analyzed 
by Booker and Bretherton (1967) is one of the cornerstone of dynamical meteorology. Its 
extension to the rotating case, even in the simple set-up with constant shear, constant 
stratification and without boundaries has been therefore studied by many investigators. 
Nevertheless, as the mathematics are involved, its exact solution has never been given 
adequately over the full domain. This paper gives this solution and interprets it, partly 
using a WKB approximation. 

In the general case, the structure of the solution is as follows. A monochromatic wave 
incident at the lower inertial level (there is no reflected wave) is amplified or attenuated 
by the valve effect, depending on the Richardson number and on the orientation of the 
wavevector. Between the inertial levels the disturbance consists of two independent solu¬ 
tions (12.4611 : in the WKB as in the QG approximation, these two solutions are necessary 
to explain a non-zero EP flux. The amplihcation (attenuation) at the lower inertial level 
is compensated by an inverse attenuation (amplification) at the upper inertial level, so 
that the overall attenuation between the two inertial levels is controlled by the decrease 
of the solution across the region in a manner akin to quantum tunnelling. The analogy 
is not complete, however, in particular because there is no reflected gravity wave. 

The situation is particularly transparent for large Richardson number J ^ 1, when 
a WKB approximation provides the solution in terms of elementary functions. In this 
case, the solution between the inertial levels is dominated by an exponentially decaying 
part, analogous to that predicted in the QG approximation. An exponentially small, 
exponentially growing solution is however present and ensures the constancy of the EP 
flux; this solution can be derived by careful consideration of the Stokes phenomenon (see 
LPV12 for details). 

When J is smaller, the decay of the balanced solution with altitude is not as pro¬ 
nounced, and the valve effect becomes significant, leading to a strong dependence of the 
wave amplitude between the inertial levels on the wavevector orientation. For J < 1, in 
particular, waves with large v, experience an amplification, with an EP flux between the 
inertial levels that is larger than that of the incident wave. This establishes a connection 
between the absorptive properties of the inertial levels and the criterion J < I for inertial 
(symmetric) instability. This connection is not apparent either in the classical transmis¬ 
sion coefficient in or in the formula for the valve effect in dlj). Practically, this 

demonstrates how small non-symmetric disturbances can become very large in inertially 
unstable flow, providing there is a small external excitation. 
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Appendix A. Derivation details 

A.l. Exact Derivation of R and T 

Following LPV12, the W functions in (l2.4aD - (l2.4cl) can be expressed in terms of hyper¬ 
geometric functions F. Using Eqs. (15.5.3), (15.5.4), (15.5.6) and (15.5.7) in Abramowitz 
and Stegun (1964) (hereafter AS), we write 

= (1 + b'; c'; T"), (Ala) 

’ = (1 + 6-*"A(a, b; c; e'), = (1 + F{a", 6"; c"; C'), (A 1&) 

= (1^1 - 5'"; c"'; l^l'^), 1^® = (|^| - I)"*" |^|-^'F(a', 6'; c'; l^r^) 

(Ale) 

In these expressions, the coefficients a, b and c are given by 


1 v — fl 
a = ~T ~ —I 

4 2 


1 

4 


— i- 


IJ- 



(A2) 


and ^ = •\/j(l + v'^) — 1/4. The other coefficients are linear combinations of a, 6 and c: 
a' = 6, 5' = 6 — c + 1, c' = 6 — a -b 1, a" = a — c -b 1, = 6 — c + 1, c" = 2 — c, a”' = a, 

b'" = a — c + 1 and d" = a — 6 + 1. Then to evaluate the complex constants A, B, C and 
D we match (j2.4a|) and (I2.46P in the vicinity of ^ = 1 and (12.461) and p.4cll in the vicinity 
of ^ = — 1 along the contour displayed in Fig. [T] (see also Eq. (E3). Using (lA 1 all - (lA Ic^ 
and the transformation formula (15.3.6) in AS, we obtain 

lU~a'(^-l)*‘^-b2-*‘^^' asC^l+, (A3a) 


W^{aA + a''B) (1 - + 2-^’'{^A + /3”B) as ^ ^ 1“, 

W - 2^‘'(Aa - Ba") + (1 -b {A/3 - Bp") ^ -1+, 

W ~ T^iCa"' + Da') + (|C| - l)~"'' {CP'" + Dp') f ^ -1". 
In these expressions the a’s and the P’s are related to the o, b, c’s by 
r(c)r(a-b 6 — c) r(c)r(c—a —6) 

- and p = 


F(a)r(6) 


F(c — a)r(c — b) ’ 


(A 36) 
(A 3c) 
(A 3d) 

(A 4) 


where F is the gamma function (see AS, chapter 6). As a first step to determine A, B, C 
and D it should be noticed that the 8 coefficients a’s and P’s can ultimately be expressed 
in terms of the 4 coefficients a, a', P', and a", given explicitly as 


, ^ r(l - ifi)r{-iv) ^ F(1 - in)T{w) 

' r(6)F(l-6*) ’ ^ r(l-a)F(a*)’ 


(A 5a) 


a = 


—2^/^^T{—iv) 


a = 


3v^r(—w) 


(A 56) 


r(a)r(6) ’ 4F(1 - a*)r(l - 6*)’ 

through the relations P = a*, P" = a"*, a"' = P'* and P"' = a'*. 

Requiring that the continuations match (I A 361) and (lA 3dl) determines the constants 
A, B and C, D through Matching near ^ = 1 and ^ = — 1 yields the linear systems 
(12.81) - (12.91) . Solving (12.811 for A and B gives 


A = 


a''P' a''*a' 


Si Si 


B = ^ 

di 6i ■ 


(A 6a) 
(A 66) 
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Solving (I2.9|) for C and D gives 

^ f'{Aa - Ba") a’{Aa* - Ba''*) 

O — --6 , 
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(A 7a) 


^ f'*{Aa* - Ba"*) ^_^^ a'*{Aa - Ba") 

bi b2 

In these equations, di and b^ are determinants appearing in explicitly given 

by 

(5i = d* a' — ad'* and <52 = \P''^ — \d\'^. (A 8) 

A closer examination of (IXSll^dXTl) shows that the Gamma functions involved appear 
through products expressible in terms of ordinary trigonometric functions using the re¬ 
flection formula 


r(z)r(i-z) = - 


SmTTZ 


(A 9) 


and the two related formulas (6.1.29) and (6.1.31) in AS. We now sketch the necessary 
calculation. For brevity we introduce the notation s = sin(7ra*) and t = sin(7r6). Apply¬ 
ing (IA9p with z = a* and z = b thus reduces the right-hand sides to 7rs“^ and 
respectively. 

Substituting (IA5I) into the two equations in (jA8l) gives 

W'-I'l' (A 10) 


Ji = 


3 ts* - st* 


and b 2 = 




2 1 / sinh(7rz/) ^ i/ sinh(7rz/) sinh(7r/i) 

These equations can be further simplified using common trigonometric identities and the 
definitions of a, b and c in terms of v and /x. A direct calculation gives 

ts* — st* = i sinh{TTv) and |s|^ — |t|^ = — sinh(7r/x) sinh(7rz/), (All) 

reducing the determinants to 


bi = —- and b 2 = — — ■ 
2v V 


(A 12) 


We can use this approach to obtain exact formulas for the transmission and reflection 
coefficients. These coefficients are defined by T = D~^ and R = CD~^, which requires 
the calculation of C and D. Substituting (IA6aP and (IA6&I1 into (I A 7611 gives 

2Re{aa"*){\l3'\^e-^'' + \d\^e^'') - ARe{aa"p'a'*) 


D = 


bib2 


(A 13) 


where the functions involving the constants d, /3', a and a" in the numerator can be 
expressed in terms of s and t as 

st* 

^ n c- 


st* - ts* 


ad'P'd* = -bib 2 




(|s|2-|t|2)(ts*-st*)’ 


(A 14a) 
(A 146) 

\s\2_\t\2 — I --"|s|2_ |i|2- (A 14c) 

Substituting (lA 141) into (lA 131) . we note that the factor ^ 1^2 appears in both the numer¬ 
ator and the denominator. Using trigonometric identities, we finally obtain the simple 
form 

D = 


= b2 




and |/3T=^2| 


(A 15) 
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or equivalently, T = D~^ = 

The derivation of a formula for the reflection coefficient follows the same argument. 
An expression for C is obtained by substituting l(X6t into (lA 7a|) . and the product of 
coefficients a', P', a" is expressed in terms of the trigonometric functions s and t, using 
the reflection formula. Hence, we deduce 


C = 26^^62^ cosh( 7 rj^)Re(a*a")] ■ 

along with the two relations 

//,o/2, /2/ //N* StT 2r2(l-fM) |sp + |tp 

aa p + a [aa j = —— 


P{a,b) {ts* — st*y 


o 2tt 2 ^^{1-ip) st*+s*t 

2a ;? Re(o 0 ) = f)» ' 


(A 16) 

(A 17a) 
(A 176) 


where P{a,b) = r(o*)r(6)r(l — 6*)r(l — a). After some direct algebra, we find 

C = 0, (A 18) 

or equivalently, i? = 0, which is valid to all orders in p. 

A.2. EP flux jump between the inertial levels 
From (j2.46ll . (I A 161) . and the definition of the Gauss hypergeometric series (15.1.1 in AS), 
we can show that near ^ = 0 

■2 1\ 

3 2 24 / 


VF = A ( 1 - - ( ^ + y ) ^2 _ -J, 


+ 7T7kM+^^' + 0(C")- (A 19) 


From evaluated very near ^ = 0 using (I2.11|) . we deduce that F^ = = 

3i {BA* — AB*) /2 between the inertial levels. Similarly, using the asymptotics form for 
IT as ^ —oo, we deduce that = —pDD* below ^ = —1. The ratio of the EP flux 

between the inertial levels to the incident EP flux is therefore given by 

eAii) 3i BA* - AB* 

Im* ■ ^ 

According to (lA 15)1 . the denominator of (lA 201) can be simplified to DD* = , while 

the numerator takes a more complicated form. Using (IA6II . we obtain 

1 


BA* - AB* = 


l<5i| 


(e2--|a'|2 - |/3f) (aa"* - a*a"). 


Upon substituting l|A 14al) and (I A 14 el) into (lA211) . we obtain 

BA* - AB* = (s*t-st*)(e2-|t|2-|,p) ^ 

3 sinh^(7ri2) sinh(7r/x) 

This equation can be further simplified using (lA 11|) and the identity 


^277Iv- U|2 


1 


|t|2 _ |5|2 = sinh(27rz/). 


Using this, (IA20I) reduces to 


^ cosh(7rt/) ^ e2^^ + 1 

|F^(i)| sinh(7r/i) e2'^^ — 1 


(A 21) 


(A 22) 


(A 23) 


(A 24) 
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A.3. Continuation of the WKB solution 

At the order fi, substitution of the WKB solutions (13.31) into expression (I3.1bl) for W{f) 
when f > 1 gives 

where cosh~^^ = In Xo continue this solution below ^ = 1, we write 

(^— 1) = (1 —^)e“®^ so cosh“^ ^ becomes In — i^l — = zsin“^ ^ —i7r/2. Therefore 

the continuation of reads 

_ ^^-l/4+w/2gi/7i-/2+i7i-/4^^ ^^-l/4-j!v/2g+/i7r/2-/ism-i (A 26) 

To continue this function below ^ = —1, we write ^ +1 = |^ + l|e“*'^ (and also ^ = —I'd)! 
so sin“^ ^ transforms into — iln _ x^ — 7 r/ 2 , leading to the continuation 

^ _^('X _ ^)-l/4+ii^/2g-i77/2|^ ^ (A 27) 

The solutions in (I3.4al) - (l3.4cl) are just (lA 25I) - (IA 271) multiplied by ie~'^^ so that the 
incident wave has a unit amplitude. 
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